Derealization of a disordered bosonic system 
by repulsive interactions 
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Clarifying the interplay of interactions and disorder is fundamental 
to the understanding of many quantum systems, including superfluid 
helium in porous media 1 , granular and thin-film superconductors^, 
and light propagating in disordered media fr 8 . One central aspect for 
bosonic systems 9 ' 11 is the competition between disorder, which tends 
to localize particles, and weak repulsive interactions, which instead 
have a delocalizing effect. Since the required degree of independent 
control of the disorder and of the interactions is not easily achievable 
in most available physical systems, a systematic experimental inves- 
tigation of this competition has so far not been possible. Here we 
employ a degenerate Bose gas with tunable repulsive interactions in 
a quasi-periodic lattice potential to study this interplay in detail. We 
characterize the entire derealization crossover through the study of 
the average local shape of the wavefunction, the spatial correlations, 
and the phase coherence. Three different regimes are identified and 
compared with theoretical expectations^*^ an exponentially local- 
ized Anderson glass, the formation of locally coherent fragments, as 
well as a coherent, extended state. Our results illuminate the role of 
weak repulsive interactions on disordered bosonic systems and show 
that the system and the techniques we employ are promising for fur- 
ther investigations of disorde red s ystems with interactions, also in 
the strongly correlated regimd^^l 

The interplay of disorder and interactions lays at the heart of the be- 
haviour of many physical systems. Notable examples are the transitions 
to insulators observed in superconductors and metals 2 5 , quantum Hall 
physics 22 , electrical conduction in DNA 23 , or light propagation in non- 
linear disordered medial An important step towards their full com- 
prehension is understanding disordered bosonic systems at zero temper- 
ature, where a competition between disorder and weak repulsive interac- 
tions is expected. Indeed, while disorder tends to localize non-interacting 
particles giving rise to Anderson localization 24 , weak repulsive interac- 
tions can counteract this localization in order to minimize the energy. 
Eventually, interactions can screen the disorder and bring the system to- 
wards a coherent, extended ground state, i.e. a Bose-Einstein condensate 
(BEC). In many years of research, mainly theoretical predictions have 
been made about the properties of the complex phases expected to appear 
as a result of this competition 9 !^!. A systematic experimental study has 
so far not been possible, since on the one hand interactions in condensed 
matter systems are strong but difficult to control 1 , while on the other 
hand in photonic systems only non-linearities corresponding to attrac- 
tive interactions2Mlhave been explored in experiments. Instead, ultracold 
atoms in disordered optical potentials are a promising system for such 
investigation*^ 2 ^! and have already enabled the observation of Ander- 
son localization for bosons in the regime of negligible interactions^! 2 !! 
Using one of these systems in a disordered lattice, we characterize the 
whole crossover from the regime of disorder-induced localization to that 
of Bose-Einstein condensation by tuning repulsive interactions in a con- 



-int 



E,„, „.*' T= 



Figure 1 | Derealization in a quasi-periodic potential. Cartoon of the 
interaction-induced derealization, a, In a very weakly interacting system with 
sufficiently large disorder, the eigenstates are exponentially localized, and sev- 
eral of the lowest energy states are populated an average of 4.4 lattice sites 
apart (Anderson glass), b, The energies of different states can become degen- 
erate due to repulsive interactions and their shape might be modified, giving 
rise to the formation of locally coherent fragments (fragmented BEC), though 
global phase coherence is not restored until c, the entire system forms a co- 
herent, extended state (BEC) at large interaction strengths. 



trolled manner. The simultaneous measurement of localization proper- 
ties, spatial correlations and phase coherence properties, and the com- 
parison with the predictions of a theoretical model allow us to identify 
the different regimes of this derealization crossover. 

The system employed consists of a three-dimensional degenerate 
Bose gas of 39 K in a one-dimensional quasi-periodic potential, which 
is generated by perturbing a strong primary optical lattice of periodicity 
d = n/ki with a weak secondary lattice of incommensurate periodicity 
n/k 2 (k = 2n/A, where A is the wavelength of the light generating the 
lattice). The corresponding Hamiltonian is characterized by the site-to- 
site tunnelling energy J of the primary lattice, which is kept fixed in the 
experiment, and the disorder strength A. The interatomic interactions 
can be controlled by changing the atomic j-wave scattering length a by 
means of a Feshbach resonance 2 ^, which in turn determines the mean 
interaction energy per particle E im (see Methods). 

In the case of non-interacting atoms, such a system is a realization of 
the Aubry-Andre model 29 , which shows an Anderson-like localization 
transition for a finite value of the disorder A/ J = 2. Above the tran- 
sition, the non-interacting eigenstates of the potential are exponentially 
localized due to the quasi-periodic perturbation of the lattice on-site en- 
ergies and the energy spectrum is split into "minibands' l 13 ! 30 !. The local- 
ization properties in this case have been studied experimentally in detail 
in ref.[27] where it was seen that several low-lying eigenstates, separated 
on average by d/(J3— 1) ~ 4Ad, where /? = k%/k\, are typically populated 
in the experiment. Adding weak interactions, the different regimes that 
appear as a result of the interplay of disorder and interactions can be ex- 



1 




-2 -1 1 2 2 4 6 

Momentum (k ) Position (lattice sites) 



Figure 2 | Analyzing the momentum distribution, a, Typical momentum 
distribution and c, mean local shape of the wavefunction recovered from a 
Fourier transform (FT) in the localized regime, and b,d in the extended regime. 
The root-mean-squared width of the momentum distributions and the exponent 
extracted from a fit (red and blue lines) to the FT give the localization proper- 
ties. The coherence properties are extracted by measuring the fluctuations of 
the phase of the interference pattern in the momentum distribution, or by the 
relative height of the two states 4.4d apart, which can be related to the spatially 
averaged correlation function g(4.4d). 

plored. For very weak repulsive interactions, the occupation of several 
eigenstates in the lowest miniband is favoured (Fig. la). This regime, 
in which several exponentially localized states coexist without phase co- 
herence, is often identified with an Anderson glass" 15 (AG). As E int is 
increased, coherent fragments, that extend over more than one well of 
the quasi-periodic potential, are expected to form (Fig. lb). In this case, 
global phase coherence would not yet be restored, and the local shape 
of the states might be modified. Some authors have called this regime a 
'fragmented BEC 12 (fBEC). Finally, for large enough E int a single, ex- 
tended phase-coherent state is expected to be formed, i.e. a macroscopic 
BEC (Fig. lc). 

The system is prepared by first loading an interacting condensate 
adiabatically from the ground state of a harmonic trap into the quasi- 
periodic lattice. The interaction energy is then slowly changed to its 
final value E int , while the confining potential is reduced. This process 
is adiabatic for most of the parameter range explored until E im becomes 
sufficiently low for the system to enter the fully localized regime. Here, 
several independent low-lying excited states are populated even when it 
would be energetically favourable to populate just the ground state. This 
loss of adiabaticity is seen experimentally as a transfer of energy into the 
radial direction (see Supplementary Information). 

The system is characterized in detail by analyzing its momentum 
distribution, which is recovered by taking an image of the condensate af- 
ter a long ballistic expansion without interactions (see Methods). From 
the momentum distribution and derived Fourier transforms, of which we 
show examples in Fig. 2, we extract the local shape of the wavefunction, 
spatial correlations, and phase coherence properties for different values 
of A/7 and E in ,IJ. The system can be approximately described as the su- 
perposition of states with the same envelope separated by 4.4d. First, the 
mean extension of individual states can be quantified by measuring the 
root-mean-squared width of the momentum distribution (Fig. 2a-b). A 
smaller (larger) width indicates a more extended (localized) state. Next, 
the mean local shape of the wavefunction on a length scale of 4.4d is ex- 




Figure 3 1 Probing the interaction-induced derealization, a, Root-mean- 
square width (in units of k\) of central peak of momentum distribution. The 
white line gives E M = 0.05A (standard deviation of the lowest-lying energies), 
where we expect the centre of the derealization crossover, b, Average expo- 
nent a of states occupying the potential wells. The line is the same as in a. 
The data taken at different values of A/ 7 and £,,„// are linearly interpolated; 
the color indicates the mean value of the measured quantity. 

traded from the Fourier transform of the square root of the momentum 
distribution. From a fit to a generalized exponential function, the local- 
ization exponent a is recovered (see also Methods), as shown in Fig. 2c- 
d. The measured momentum width and exponent are shown in Fig. 3. We 
find that for very small £,„,, the states are exponentially localized, since 
a « 1, and the momentum width is large, consistent with the Anderson 
glass regime. Increasing E jnt , the width decreases while the exponent 
increases up to or ss 2. Repulsive interactions therefore delocalize the 
system as expected, or alternatively, the localization transition is shifted 
to higher values of the disorder strength A/7 when interactions are intro- 
duced into the system. The position of the derealization crossover is in 
good agreement with the expectations of a simple screening argument 
The increasing interaction energy serves to smooth over the disordering 
potential in the occupied sites, providing a flatter energetic landscape 
on which more extended states can form. The centre of the crossover is 
therefore expected to occur when E- mt is comparable to the standard devi- 
ation of energies in the lowest miniband of the non-interacting spectrum, 
0.05A (white line in Fig. 3, see also Supplementary Information). 

The correlation properties of neighbouring states can be extracted 
from a Fourier transform of the momentum distribution itself, which 
gives the spatially averaged first order correlation function g(x) (see 
Methods). In Fig. 4a-b, g(x) at x = 4.4d is shown for both the experiment 
and a ground-state theory that we have developed, with generally good 
agreement. In the localized regime, the correlation is exactly zero in the 
theory, since no neighbouring states are occupied. In contrast, the corre- 
lation is finite in the experiment due to the occupation of neighbouring 
localized states arising from the non-adiabatic loading, but is small since 
the states are independent. As E im is increased, the correlation features 
a crossover towards larger values, signalling that coherence is progres- 
sively established locally over distances of at least 4.4d. The shape of 
the crossover in the experiment is again in qualitative agreement with the 
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Figure 4 | Probing the phase coherence of the system, a & b, Correlation of neighbouring localized states, a, Experiment (white line as in Fig. 3), b, 
Theory of the ground state. The white and orange dashed lines show the phase boundaries between an Anderson glass (AG), a fragmented BEC (fBEC), and 
a macroscopic Bose-Einstein condensate (BEC). c, Correlation of neighbouring states for A/7 = 12, corresponding to black dashed line in a. The error bar is 
given by the standard error of the mean, d, Standard deviation of the phase measured by repeating the experiment up to 25 times for a given set of parameters for 
A/7 = 12. The error is estimated as A0/ VA', where N is the number of images from which the phase was extracted. The blue bar shows the phase fluctuations 
measured for an extended system below the localization threshold. The grey dash-dotted line gives the standard deviation for a purely random distribution. In 
C & d, the grey dashed (dotted) lines give £,,„ = 0.05A (0.17A). 



expectation of the screening argument above (see white line in Fig. 4a). 

Finally, information about the phase coherence of neighbouring 
states can be obtained by measuring the phase <p of the interference pat- 
tern in the momentum distribution for repeated runs of the experiment 
with the same parameters (see Methods for details). If the states are not 
phase locked, <p changes almost randomly at each repetition of the ex- 
perimental sequence. In Fig. 4d we show the standard deviation of (p, 
estimated from a large number of repetitions of the experiment, for fixed 
A/7 = 12. We see a slight decrease of the phase fluctuations with increas- 
ing £,„,, that nevertheless remain relatively large in the crossover region 
where the correlation increases (Fig. 4c). The fluctuations finally drop 
to the background value only when E jm is comparable to the full width 
of the lowest miniband of the non-interacting spectrum, 0.17A. These 
observations confirm that in the localized regime the states are totally 
independent, which together with the localization properties (Fig. 3) in- 
dicates that the system can indeed be described as an Anderson glas d 11 ! 15 !. 
The system crosses a large region of only partial coherence while becom- 
ing progressively less localized as £„„ is increased. This is consistent 
with the formation of locally coherent fragments expected for a fBEC. 
An analogous fragmentation behaviour was reported in ref . [31] Finally, 
the features of a single extended, coherent state are seen, i.e. a BEC. 

In the mean-field theory, boundaries between the different phases ex- 
pected for the system can be defined (see Methods). In particular, the 
transition from the Anderson glass phase to a fragmented BEC (white 
line in Fig. 4b) occurs when g(AAd) starts to increase. Similarly, the 
orange line in Fig. 4b shows where the fragments are locked together 
in phase to form a single macroscopic condensate for very large inter- 
actions. The generally good agreement between the experimental ob- 
servables and theory indicates that our system is well described by the 
mean-field theory for most of the parameter space explored experimen- 
tally. 

In conclusion, we have provided the first experimental characteriza- 
tion of the localization, correlation and coherence properties of the vari- 
ous regimes due to the competition of disorder and weak repulsive inter- 
actions in a bosonic system. Other aspects of the derealization crossover 
worth further study are, e.g. the detailed properties of the ground state of 
the AG regime, which was not possible to study in the present set-up, 
and the presence of a superfluid-insulator transition at the BEC-fBEC 
boundary analogous to the one observed in superconductors 2 . Regard- 
ing the latter, in transport experiments analogous to the ones described 
in ref. [27] we have been able to verify that the AG and fBEC regimes 
are not inconsistent with being insulating, as is the case in the regime 



of vanishing E mt (see Supplementary Information). Finally, it would be 
appealing to employ the present system and the correlation analysis in- 
troduced here to explore the regime of strong correlations, E int » 7, 
which could be reached by using a quasi- ID system with strong radial 
confinement. There, another elusive insulating phase due to the coop- 
eration of disorder and interactions, the so called Bose-glass phase, is 
expected to a ppear, alth ough there is debate on the exact shape of the 
phase diagra m 11 ! 15 ! 18 ! 19 !. 

METHODS 

Condensate with tunable interactions. A K condensate of about N = 
20, 000 atoms with an .s-wave scattering length of 250ao, where ao = 52.9 pm 
is the Bohr radius, is prepared in a harmonic optical trap. The condensate 
is loaded into the quasi-periodic potential while the optical trap is decom- 
pressed in about 250 ms to reduce the harmonic confinement, and a gravity- 
compensating magnetic field gradient is added. At the same time, the scat- 
tering length a is changed by means of a broad Feshbach resonance to values 
ranging from a < 0.1 ao to about a = 300ao (ref. 1281 

Quasi-periodic potential. The quasi-periodic potential is created by two ver- 
tically oriented laser beams in standing-wave configuration. The primary lat- 
tice is generated by a Nd:YAG laser with a wavelength of A\ = 1064.4 nm 
and has a strength of si - Vi/£rj = 10.5 (corresponding to J/h = 79 Hz), 
as measured in units of the recoil energy Eg_ i = h 2 /(2MA 2 ). The secondary 
lattice is generated by a Ti:Sapphire laser of wavelength A2 = 866.6 nm, the 
strength being adjustable up to S2 = Vi/Eri = 1-7. Both beams are focussed 
onto the condensate with a beam waist of about 150 pm. The lattice lasers give 
a harmonic confinement of i±> ± = 2n x 50 Hz in the radial direction. In the ver- 
tical (axial) direction, a weak confinement of 5 Hz is given by a weak optical 
trap as well as by a curvature from the gravity-compensating magnetic field. 
Energy scales. In the tight-binding limit, the hopping energy 7 and disor- 
der strength A can be estimated as 7 = 1.43.?° 98 exp J-2.07 yTj"} £r i and 

A = 0.5i 2 j8 2 (l.0264exp(-2.3624/^- 59 })£: R ,i (ref.fJO). The experimental un- 
certainty on A/7 is around 15%. We estimate that around 30 lattice sites, 
corresponding to about 7 localized states, are populated during the load- 
ing of the lattice. We then define a mean interaction energy per particle 
Ej„i = gN/7 J|i/;(r)| 4 d 3 r, where g = Ajitp-a/m and ip(r) is a Gaussian ap- 
proximation to the on-site Wannier function. We include coupling into the 
radial directions of our system, with the consequence that the interaction en- 
ergy is non-linear in the scattering length. Though this definition of the energy 
is strictly valid only in the localized regime, comparison with a numerical 
simulation of our experimental procedure has shown that it is a good approxi- 
mation for all values of the scattering length up to an error of 30%. Note that 
the potential energy from the residual harmonic confinement is approximately 
3 x 10~ 3 7 over a distance 4.4d. 
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Momentum distribution analysis. The images of the momentum distribu- 
tion are taken by absorption imaging with a CCD camera after 36.5 ms bal- 
listic expansion. At the time of release, the scattering length is set to below 
lao in less than one ms and kept there until the Feshbach magnetic field is 
switched off 10 ms before taking the image - at this point, the system has 
expanded a sufficient amount to minimize the effect of interactions. For such 
a free expansion, the image is approximately the in-trap momentum distribu- 
tion p(k) = (T' 1 (ky¥(k)) fref. 1211 The acquired images are integrated along 
the radial direction to obtain a profile. In momentum space, the width of the 
central peak is calculated by taking the root mean square width within the first 
Brillouin zone. Due to the quasi-periodic lattice potential, for a sufficiently ho- 
mogeneous system the in-trap wavefunction can be decomposed into copies of 
a single state with real and non-negative envelope %(x) ~ exp(-|x/L| ff ), spaced 
by 4.4d. Therefore in momentum space, ^jp(k) = ^(k)S{k), where S(k) is 
an interference term, and %(x) can be extracted from a Fourier transform of 
yfpik) (see also Supplementary Information). We fit to the sum of two gener- 
alized exponential functions A exp (-|(x - x c )/L\ a )- [1 + B cos(ki(x- x c ) + 6)], 
where x c denotes the centers of each of these functions, spaced by 4.4a!. 
From this fit, the exponent a is recovered. In addition, from the Wiener- 
Khinchin theorem, the momentum distribution can be expressed in terms 
of the first order correlation function G(x',x + x') = (1T(x?y¥(x + x')) as 
p(k) oc g" 1 jG(x', x + x')dx' . By taking the Fourier transform of the momen- 
tum distribution itself, we can therefore recover the spatially averaged corre- 
lation function g(x) = j G(x' ,x + x')dx' . With the same fitting function as 
above, we evaluate the spatially averaged correlation between two states 4.4 
lattice sites apart, Aj/Ai. Experimentally, the correlation function saturates at 
a value around 0.5 due to the finite momentum resolution. The fluctuations 
in phase between neighbouring states are seen as a fluctuation of the phase 
(f> of the interference pattern of the momentum distribution, which is directly 
extracted from a fit (see also Supplementary Information). The 2D graphs in 
Figs 3 and 4 were generated by linearly interpolating a total of 1 30 averaged 
datapoints at 9 different values of disorder, changing the interactions. Typical 
experimental scatter and statistical errors are seen in Fig. 4c. 
Theory of the ground state. The theoretical calculations presented in the 
paper rely on a meanfield approach similar to the one of ref. [32] This is an ef- 
fective one-dimensional model which partially includes also the radial to axial 
coupling, and is known to provide an accurate description in the two limiting 
cases of Anderson localization and BEC. The boundaries between the different 
regimes shown in Fig. 4b are obtained by analyzing the correlation function 
g(x) and the density distribution. In the theory we define the Anderson glass 
phase as the one in which the correlation g(4.4d) is zero. To enter the frag- 
mented BEC (fBEC) phase, we require g(4.4d) > 0, which implies that co- 
herent fragments composed of adjacent localized states can start to form. For 
increasing E mt the extension of the fragments increases, until most of the sys- 
tem remains in a single component, which corresponds to a macroscopic BEC. 
To define the boundary between the fBEC and the BEC regimes, we first iden- 
tify as fragments the parts of the system separated by low-density regions for 
which an applied relative phase twist does not affect the energy of the system. 
When one single macroscopic fragment forms, we assume the system to be in 
the BEC regime. A more detailed description of the theoretical methods can 
be found in the Supplementary Information. 
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Supplementary Information 



Energy spectrum The potential of our system can be written in general 
as 



V(x) = s 1 E f!i sin 2 (kix) + S2E R2 sin 2 (k 2 x + </>) + V ext (x, r ± ), 



(SI) 



where E RJ = % 2 k 2 j{2M) = h 2 /(2MA 2 ) is the recoil energy for the lattice 
with wavelength /i, = 2n/kj, and s f = V,/£r_, is the height of the lattice 
( in units of £r_, . Each of the two lattices is A,/2 periodic. Any external 
confining potential is given by V exl (x,r ± ). The lattice spacing of such a 
potential is to good approximation d = X\j2. If the ratio /? = k 2 jk\ is 
an irrational number, eq. lSll describes a quasi-periodic potential. In our 
case, X\ = 1064.4 nm and A 2 = 866.6 nm, giving » 1.228. 

The essential features of such a potential are visible in Fig. IS 1 1 The 
potential energy minima of the primary lattice are modulated by the 
second one, giving rise to characteristic wells separated on average by 
1/08 — 1) ~ 4.4 lattice sites. The energy scales that characterize the cor- 
responding Hamiltonian to such a potential are the tunnelling energy of 
the primary latticed 



7 = 1.43s?- 98 exp {-2.07 V*T}£r,i. 
and the disorder energ> 2 

A = 0.5s 2 f3 2 (1.0264 exp {-2.3624/ j?- 59 }) £ R ,, 



(S2) 



(S3) 



Neglecting the external confining potential, the spectrum of such a 
quasi-periodic potential can easily be calculated and is shown in Fig. IS2I 
for various values of the disorder strength A/ 7. A striking feature is the 
appearance of minigaps in the spectrum, the lowest of which has approxi- 
mately the same width for all values of A/7. A minigap appears when the 
potential has two neighbouring lattice sites with almost the same mini- 
mum potential energy. Locally, the potential then looks like a double 
well, for which the two lowest-lying eigenstates have an energy split- 
ting of 27. In fact, the width of the lowest minigap is approximately 2J 
throughout the range of A/7 shown. The lowest "miniband" of energies 
corresponds to the lowest energy eigenstates localized in the potential 
wells 4Ad apart. Since in the experiment, only the states in the first 
"miniband" are populated, we restrict our analysis to these energies and 
find that their standard deviation is approximately 0.05A, while the ex- 
tension of this band is approximately 0.17A. The effect of a confining 
potential on the spectrum has been analysed previously in ref.[2] 




Position 

Figure S1 | Quasi-periodic potential. The quasi-periodic potential realized 
in the experiment for lattice incommensurability /? = 1.228 . . . and A/7 = 15. 
The red stripe shows the amplitude 2A of the perturbation by the secondary 
lattice. The blue stripe denotes the position of the first Bloch band of the 
primary lattice, with width 47. The quasi-periodic lattice is characterized by 
potential wells approximately every 1/(J3 - 1) ss 4.4 lattice sites, which arise 
from the beating of the two lattices (dashed grey line). 




Figure S2 | Energy spectrum of the quasi-periodic potential. Energies 
of eigenstates of the non-interacting system as a function of A/ 7. 



Momentum distribution analysis Fourier transform techniques are 
used to extract information both about the local shape of the wavefunc- 
tion, and about the coherence properties of neighbouring states. Af- 
ter a long free expansion without interactions, the image of the atoms 
that is acquired is approximately the in-trap momentum distribution 
p(k) = OF* (&)¥(£)), where is the Fourier transform of the bosonic 
field operator ^(x). In order to recover information about the in-trap 
wavefunction, we can therefore use an inverse Fourier transform. 

Due to the quasi-periodic nature of the employed lattice potential, 
we expect that for a sufficiently homogeneous system, the in-trap wave- 
function can be decomposed into copies of the same state with real and 
non-negative envelope £(x), spaced by D = 4Ad. The overall wavefunc- 
tion can therefore be approximated as 



Mx) = ^] a^(x - jD)e 



-•tpj 



(S4) 



where 0^ is the local phase, and a useful example of £(x) is a generalized 
exponential function exp(-|x/L|"). In momentum space, the magnitude 
of the overall wavefunction can then be written as -Jp(k) = \g(k)\S(k), 
where 



S(k). 



-i(jkD+<Pj) 



(S5) 



is an interference term. For many envelope functions £(x), such as the 
generalized exponentials with < a < 2, the Fourier transform ^(k) it- 
self is real and non-negative 3 , so that the inverse Fourier transform of 
^jp(k) can be written as £(x) o S(x). This is simply the convolution of 
the envelope of a single state f (x) with the Fourier transform of the inter- 
ference term, <S(x), which can be approximately described as a series of 
sharp peaks (approaching iS-distributions) spaced by D, with a decreasing 
amplitude and phases that depend on the local phases <j>j and amplitudes 

The inverse Fourier transform of the square root of the momen- 
tum distribution p{k) therefore gives the average local shape of the 
(wave)function £(x). Due to our finite resolution in momentum space 
(about k\ 120), we are only able to resolve easily two neighbouring states. 
The averaged wavefunction is analysed by fitting to the sum of two gen- 
eralized exponential functions modulated by the primary lattice 



Yj A/ exp (-|(x - xj)/L\") [l + B cos(£i (x - xj) + Sj)] 



(S6) 
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where X\ = and x 2 = A Ad (see Fig. [S3] for examples). From such 
a fit, the exponent a can be extracted. The local size L is shown in 



5 



aWa ' \ 




a 


" *. *.* 


\ " v 


b 



1 2 3 4 5 6 7 

Position (lattice sites) 



Figure S3 | Examples of Fourier analysis. Typical Fourier transforms 
(dots) of the square root of the momentum distribution giving the average lo- 
cal shape of the wavefunction, a, in the localized regime, b, for a coherent, 
extended state. The red and blue lines give the fits to the data for a generalized 
exponential function centered at the origin (red line) and at 4.4c! (blue line). 



Fig. IS4I which shows a complementary behaviour to the momentum 
width (Fig. 3a in the main text). 

On the other hand, the inverse Fourier transform of the momentum 
distribution itself can be employed to find the correlation properties of 
neighbouring states. Using the Wiener- Khinchin theorem, the momen- 
tum distribution p(k) can be related to the first order correlation func- 
tion G(x\x + x') = <$*{x , y${x + x')>, as p(k) oc g" 1 / G(x',x + x')dx'. 
By taking the Fourier transform of the momentum distribution, we can 
therefore recover the spatially averaged correlation function g(x) = 
jG(x',x + x')dx'. We fit with the same generalized exponential of 
Eq.[S6]and recover the spatially averaged correlation between two states 
4.4 lattice sites apart as A 2 /A[. Also here, the finite momentum resolu- 
tion limits our analysis to two neighbouring sites, and it follows that the 
correlation function saturates at a value around 0.5. 

The effect of a fluctuating phase between neighbouring states is seen 
as a shift of the phase <p of the interference in the momentum distribution. 
We extract this phase by fitting the momentum distribution directly with 
a fitting function 

^ A- k txp{-{k-~kf 1 2w 2 )-[\+B cos [D(k-~k) + (ST) 

where k c is the center of the distribution, determined by fitting the aver- 
age of all images of a given dataset. 

Adiabaticity As discussed in the paper, in the experiment the system is 
prepared by first loading a strongly interacting condensate (a = 250ao) 
into the quasi-periodic potential, in presence of a tight axial harmonic 
confinement (60 Hz). The scattering length is then reduced to its final 
value with an exponential ramp lasting 250 ms (r as 25 ms), while the 
harmonic confinement is linearly reduced to a much lower value (5 Hz). 
We check the degree of adiabaticity of this procedure by monitoring the 
evolution of the radial degrees of freedom. From the radial profiles ex- 
tracted from the absorption images, we can measure the radial temper- 
ature of the system. We find that there is always either a quasi-pure 
condensate or condensed and thermal components. In the latter case, 
the temperature can be directly measured from the thermal component, 
while in the former case, only an upper bound for the temperature (about 




Figure S4 | Local size of States. The local size L, as extracted from fits to the 
local shape of the wavefunction, gives the average extension of the individual 
states. This quantity is strictly valid only in the localized regime. The white 
line gives the expected position of the delocalization crossover (£,„, = 0.05A). 




Figure S5 | Adiabaticity of preparation. The radial temperature is esti- 
mated from radial profiles extracted from the images. The white line gives 
the expected position of the delocalization crossover (£,„, = 0.05A). The 
black crosses show the positions in the A/J-Ei,„/J plane of the datapoints 
from which all 2D graphs were interpolated. 



20 nK) can be given. The condensation temperature T c varies across 
the parameter range explored, but we estimate^ that it is on the order of 
60-100 nK. 

We interpret the radial excitation as an interaction-mediated transfer 
of excitation energy from the axial to the radial degrees of freedom. What 
presumably happens in the system is the following: As the tunnelling 
time between the states separated by 4.4d becomes longer than the exper- 
imental timescale for decreasing the interaction strength, the population 
of these states cannot move adiabatically to the absolute ground state of 
the system. The states are therefore left with an excess population, hence 
an excess interaction energy that spatially broadens them by populating 
even more excited states in the neighbouring sites of the primary lattice. 
These excited states have an energy on the order of the mean separation 
of energies of neighbouring lattice sites 0.8A, and can decay to the low- 
est state in the potential well on a timescale of A -1 by transferring their 
excess energy to the radial directions. 

The observation that this process takes place only near the onset of 
the Anderson glass regime is in agreement with the expectation of sup- 
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pressed tunnelling between the states 4.4d apart in the exponentially lo- 
calized regime. A radial excitation therefore shows that there are axial 
excitations, e.g. due to a lack of adiabaticity during the loading proce- 
dure. On the other hand, a lack of radial excitations is not strictly a proof 
that there are no axial excitations, since the radial heating is caused by 
local axial excitations, i.e. arising from neighbouring lattice sites, while 
there could be axial excitations over longer distances that are not able to 
release energy by relaxing into the true ground state. 
Ground state phase diagram For the calculation of the theoretical pro- 
files and ground state phase diagram presented in this paper we have 
used the following meanfield approach. In the limiting cases of a sin- 
gle Anderson localized state (AL) or of a coherent condensate (BEC), 
the system can be described by a single wavef unction if/(x, r ± ) that is a 
solution of the Gross-Pitaevskii equation in the BEC regime, or of the 
Schrodinger equation in the AL regime. In the latter case the prob- 
lem is separable, and the wavefunction can be factorized as \j/(x, r ± ) = 
<p(x)(p(r ± ). In the intermediate regime of a glassy phase made of sev- 
eral independent fragments, the atomic distribution that minimizes the 
energy can be accounted for by considering a wavefunction of the form 
i//(x,r ± ) x Y*i<Pi( x )<Pi( r ±) (similarly to the approach used by Lugan et 
a/PJ). An approximate way to describe the overall behaviour of the sys- 
tem in the crossover between different regimes is to consider an effective 
wavefunction tf/(x, r ± ) = <p(x)if>(r±, cr(x)), where the radial component is 
a Gaussian with an x-dependent width, and minimizing the correspond- 
ing energy functional 



E[(p, o~] = J dxif*(x) 



T, 2 



ti 2 1 



V j + V(x) + - + 



2 2ner 2 



2in 

M 4 



2m cr 2 



ma) ± 
~2~ 



ip(x) 



(S8) 



with V(x) given by eq. IS 1 1 that corresponds to solving the non- 
polynomial Schrodinger equation of ref. [6] This approach has the ad- 
vantage of capturing the modification of the axial density distribution in 
the crossover from the localized to the extended regime, and we use it to 
evaluate the spatial and momentum distributions of the ground state of 
the system (Fig. lS6t . In particular, from the calculated momentum dis- 
tributions we extract the correlation function g(4Ad) shown in Fig. 4b as 
we do for the experimental data. 

The boundaries between different regimes in Fig. 4b have been de- 
rived by studying both the g(4.4d) and the spatial distributions. For ex- 
ample, the system is expected to pass from the Anderson glass (AG) to 
the fragmented BEC (fBEC) when neighbouring localized states start to 
be macroscopically occupied and have a relative phase coherence. On 
the one hand, a non-zero value of g(4.4d) indicates the occupation of 
neighbouring states. On the other hand, we can check for the coherence 
properties of such states by studying the effect of introducing a phase 
twist over an extension d between them. If the two states are indepen- 
dent, the phase twist will not modify the energy of the system, while the 
contrary happens if they are coherent. In practice, we evaluate the energy 
cost per particle, 

h 2 (2tt) 2 



6E k 



-lvte)r 



(S9) 



2m d 

for introducing a 2n phase twist at the lattice site k where the two states 
connect. If SEt is less than the change in the energy per particle associ- 
ated with the removal of one atom from the system, E jm /N, where E mt 
is the last term in eq. IS8I then the two states are considered indepen- 
dent. We define the boundary between the AG and fBEC regimes where 
g(4.4d)=0.0l. We have verified that above this threshold the neighbour- 
ing states are also phase locked, i.e. they constitute a coherent fragment. 
For increasing interaction energy the size of these coherent fragments in- 
creases, while their number decreases. Due to the harmonic confinement, 
the system tends to form a large central core surrounded by fragments in 
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Figure S6 | Theoretical density distribution. Density distribution calcu- 
lated from the ground state theory for A/7 = 10, a, on a linear scale, and b, 
on a logarithmic scale. The interaction energy increases from top to bottom 
(Ei„,/J = 0.003, 0.066, 0.26, 1 .5). The top two panels are for £,„, such that the 
system is an Anderson glass, the third panel shows a fragmented BEC, and 
the lowest panel depicts a BEC. The red lines show Ei„,/N, which is used to 
evaluate the different phases. 



the low-density tails. We assume that the system has entered the BEC 
regime when there is a single macroscopic fragment at the center of the 
trap, and the population of each of the outer fragments is less than 1% of 
that of the central one. 

Transport experiments In the experiment, in addition to studying the 
momentum distribution, we have investigated the transport properties of 
the various regimes of the system, using the same technique we already 
employed in ref. [7] The technique consists of suddenly releasing the ax- 
ial harmonic confinement, while keeping both the quasi-periodic lattice 
potential and the radial confinement, and then observing the subsequent 
diffusion (or lack thereof) of the atomic cloud in the axial direction. In 
the previous experiment we observed that for vanishing £,,„ the diffu- 
sion becomes strongly suppressed for A/ J > 2. We have now observed 
that it continues to be suppressed for A/J>2 also for the values of E int 
explored in the present work, irrespective of whether the system is in 
the AG, fBEC or BEC regimes. We can interpret the absence of diffu- 
sion in the AG (fBEC) regime as a result of the insulating nature of the 
system, which occupies (partially) localized states. In contrast, in the 
extended BEC regime, a slow subdiffusive expansion is expectec 8 , or the 
expansion might be completely suppressed for some values of the param- 
eters by self-trapping 9 . The slow expansion however would probably not 
be detectable on the one second timescale of the experiments we have 
performed so far. The absence of a different diffusion behaviour of the 
system in the AG, fBEC and BEC regimes unfortunately does not allow 
us to identify the boundary between the fBEC and BEC regimes using 
this technique, nor to unambiguously prove the insulating nature of the 
AG and fBEC regimes. Other methods therefore need to be developed 
for such a purpose. 



Supplementary References 



Gerbier, F. etal. Interference pattern and visibility of a Mott insulator. Phys. Rev. 
A 72, 053606 (2005). 

Modugno, M. Exponential localization in one-dimensional quasi-periodic optical 
lattices. New J. Phys. 11,033023(2009). doi:1 0.1 088/1 367-2630/11/3/033023. 
Giraud, B. & Peschanski, R. On positive functions with positive Fourier trans- 
forms. Acta Phys. Pol. B 37, 331-346 (2006). 



7 



4. Burger, S. etai. Quasi-2D Bose-Einstein condensation in an optical lattice. Eu- 
rophys. Lett. 57, 1-6(2002). 

5. Lugan, P. et at. Ultracold Bose gases in 1 D disorder: From Lifshits glass to 
Bose-Einstein condensate. Phys. Rev. Lett. 98, 170403 (2007). 

6. Salasnich, L, Parola, A. & Reatto, L Effective wave equations for the dynamics 
of cigar-shaped and disk-shaped Bose condensates. Phys. Rev. A 65, 043614 
(2002). 

7. Roati, G. etal. Anderson localization of a non-interacting Bose-Einstein conden- 
sate. Nature 453, 895-898 (2008). 

8. Larcher, M., Dalfovo, F. & Modugno, M. Effects of interaction on the diffusion of 
atomic matter waves in one-dimensional quasiperiodic potentials. Pfiys. Rev. A 
80, 053606 (2009). 

9. Anker, T. et al. Nonlinear self-trapping of matter waves in periodic potentials. 
Phys. Rev. Lett. 94, 020403 (2005). 



8 



